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SPECIFICATION 

Electronic Version 1.2.8 
Stylesheet Version 1 .0 

Discrete Multitone Modem 
Initialization System and Method 

Background of Invention 

[0001] Priority Data 

/f / 

f ^^[0002] This application is a continuation of co-pending United States Patent Application 
Serial No. 1 0/064,570, filed jLf/26, /o02. 

[0003] Field of Invention 

[0004] The present invention relates to the field of digital communications. More 
specifically, the present invention provides for a novel initialization system and 
method for discrete multitone modems. 

[0005] Background of Invention 

[0006] The demand for high-speed data rate services has been increasing very rapidly. 

One way to meet this demand is to expand the existing infrastructure by connecting 
all users of bandwidth-consuming applications to fiber optic networks. However, this 
solution is not financially viable. An attractive alternative to fiber optics involves the 
transfer of information at high-speed data rates over ordinary telephone networks 
using Orthogonal Frequency Division Multiplexing (OFDM). When applied in a wireless 
environment, it is generally referred to as OFDM, while in a wired environment the 
term Discrete Multitone (DMT) is more appropriate. 

[0007] 

The history of OFDM dates back to 1 966, when Robert W. Chang published his 
paper on the synthesis of bandlimited signals for multichannel transmission without 
interchannel interference (ICI) and intersymbol interference (ISI). Chang was awarded a 
patent for his work in 1 970. In 1 971 , S.B. Weistein and P.M. Ebert contributed to 
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OFDM with the use of Discrete Fourier Transform (DFT) to perform baseband 
modulation and demodulation, removing the need to use banks of subcarrier 
oscillators. But, Weistein and Ebert's system could not guarantee orthogonality. In 
1980, A. Peled and A. Ruiz solved the orthogonality problem with the introduction of 
the cyclic prefix (CP). OFDM is currently being used as the European Digital Audio 
Broadcasting (DAB) standard and under the name DMT, has been chosen as the 
modulation scheme for ADSL, a technique for transmission over twisted pair 
telephone lines standard and is a major contender in the ongoing VDSL standard for 
high-speed transmission on the existing telephone network. 

[0008] The principles of multicarrier modulation are described more completely in " 

Multicarrier Modulation for Data Transmission: An Idea Whose Time Has Come \ by 
John A.C. Bingham, IEEE Communications Magazine, Vol. 28, No. 5, pages 5-14, May 
1990. As is known in a modem system utilizing multicarrier modulation such as DMT, 
the bits of input data for transmission within each block or symbol period are 
allocated to sub-carriers. 

[0009] At the receiver DMT signals are subject to synchronization errors due to oscillator 
impairments and sample clock differences. The effective implementation of DMT 
requires that the received signal be recovered with a near perfect synchronization. 
Two types of synchronizer structures are well known in the art. Data-aided 
synchronizers use the receiver's decisions or a training sequence to compute the 
timing offsets. Non-data-aided synchronizers operate independent of the transmitted 
information sequence. For the NDA approach, the timing estimates maximize the NDA 
likelihood function, which is obtained by averaging the likelihood function over the 
random information variables. A thorough explanation of DMT synchronization can be 
found in, " Synchronization with DMT Modulation" \ Thierry Pollet and Miguel Peeters, 
IEEE Communications Magazine, April 1 999, and in, " ML Estimation of Time and 
Frequency Offset in OFDM Systems % \ Jan-Jaap van de Beek, IEEE Transactions on 
Signal Processing, Vol. 45, No. 7, July 1997, both herein incorporated by reference. 

[0010] 

It is also necessary for the receiver to perform channel estimation, which allows 
for compensation of the effects of the transmission medium (ex. telephone line 
channel). Prior art systems for channel estimation are well known in the art. Channel 
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estimation in OFDM systems is usually performed with the aid known pilot symbols. A 
detailed explanation of channel estimation and pilot tones can be found in, " Pilot 
Tone Selection for Channel Estimation in a Mobile OFDM System ", Rohit Negi and John 
Cioffi, IEEE Transactions on Consumer Electronics, Vol. 44, No. 3, August 1998, 
incorporated herein by reference. 

[001 1] While many techniques for synchronization and channel estimation are known in 
the art, the prior art systems require complete synchronization prior to channel 
estimation. A system that efficiently combines synchronization and channel estimation 
will serve to significantly reduce the start-up time of any communication systems. 

Summary of Invention 

[001 2] The longstanding but heretofore unfulfilled need for an improved link 

initialization method for system level design of a modem using Discrete Multitone 
(DMT) technique, which enables very high bit-rate transmission and long distance 
reach over the existing telephone lines is now met by a new, useful, and nonobvious 
invention. 

[001 3] The present invention can be implemented into a modem system, and into a 

method of operating the same to initialize the operating conditions of the modem-to- 
modem session. 

[0014] The present invention discloses an improved initialization method for modem 
communication wherein the timing offset required for synchronization is estimated 
utilizing an entire received DMT frame and the channel impulse response is estimated 
utilizing the pilot tones received within the same DMT frames utilized for 
synchronization. 

[001 5] In a preferred embodiment of the invention, the timing offset estimation and the 
channel impulse response are performed simultaneously. 

[001 6] The present invention provides for the incorporation of a plurality of DMT frames 
and a plurality of pilot tones within the frames. 



[0017] 



In a preferred embodiment of the initialization method of the present invention, a 
plurality of DMT frames of length N+LP are transmitted. The DMT symbol to be 
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transmitted comprises N samples and the associated cyclic prefix has a length of LP 
samples. Multiple DMT frames are received at a receiver. Interpolation is performed on 
the received DMT frame. The integer timing offset and the fractional timing offset are 
then estimated using non-data aided maximum likelihood correlation with a pre- 
stored frame of length N+LP. For the channel estimation, a Minimum Mean Square 
Error criterion of the pilot tones received within the plurality of DMT frames is used, to 
estimate the channel impulse response. 

[001 8] The present invention also provides for correction of the integer timing offset in 
the time-domain and fractional timing offset in the frequency-domain. In one aspect 
of the invention, symbol timing offset correction is performed in the time-domain 
responsive to the integer timing offset by adjusting the symbol clock. In another 
aspect of the invention, sample timing offset correction is performed in the 
frequency-domain responsive to the fractional timing offset through the delay-rotor 
property. 

[001 9] The present invention also provides" for the synthesis of a 1 -tap frequency domain 
equalizer based on the estimation of the channel impulse response. The channel 
impulse response is padded with zeroes to accommodate for circular convolution 
prior to the synthesis of the 1 -tap frequency domain equalizer. 

[0020] In a preferred embodiment of the present invention, known symbols are used for 
the pilot tones. It is well known that DMT modulation schemes exhibit high peak-to- 
average power ratio (PAR). It is arguably one of the greatest drawbacks on DMT. High 
PAR values demand expensive and less efficient linear amplifiers at the analog front- 
end; therefore, overall system complexity and cost increases significantly. 
Constellation shaping can be used to reduce the peak power while maintaining the 
same bit error rate (BER) performance. It has been shown that through constellation 
shaping, peak power of a DMT signal can be reduced by 3-6 dB without sacrificing 
bandwidth. Optimal peak power reduction can be achieved through the use of a 
symmetric spherical QAM constellation as disclosed in co-pending application, U.S. 
Serial Number 1 0/064,570 filed on July 26, 2002. 

[0021] In another preferred embodiment of the present invention, an optimal set of 
known pilot tones are transmitted to the receiver. 
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[0022] It is an object of the present invention to provide a method and system for data 
communications in which transmitting and receiving modems may initialize 
communications in order to maximize the available data communication rates and 
accuracy. 

[0023] It is a further object of the present invention to provide for a fast link initialization 
between modems without significantly increasing complexity in the system 
architecture. 

[0024] It is also an object of the present invention to provide a method and system, 

which improves the over-all system performance through reduction in symbol errors. 

[0025] It is another object of the present invention to provide such a method and system 
that is particularly suited for modems that may be realized in a small number of 
integrated circuit devices. 

[0026] The invention accordingly comprises the features of construction, combination of 
elements, and arrangement of parts that will be exemplified in the description set 
forth hereinafter and the scope of the invention will be indicated in the claims. 

[0027] As previously discussed in the background of the invention, OFDM when applied 
in a wired environment in appropriately referred to as Discrete Multitone (DMT). While 
the detailed description of the invention is focused on the method and apparatus as 
applied to a DMT environment, it is within the scope of the invention to apply the 
methods described in an environment employing various orthogonal multi-carrier 
signaling systems to include OFDM as well as DMT. 

[0028] Other aspects and advantages of the present invention can be seen upon review of 
the figures, the detailed description, and the claims, which follow. 

Brief Description of Drawings 

[0029] For a fuller understanding of the nature and objects of the invention, reference 
should be made to the following detailed description, taken in connection with the 
accompanying drawings, in which: 

[0030] FIG. 1 is a system level diagram of the present invention. 
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[0031] FIG. 2 describes the mathematical formulae of the transmitted signals with noise 
and the channel impulse response used inpresent inventionFIG. 3 describes the 
mathematical formulae of the correlation function. 

[0032] FIG. 4 describes the mathematical calculations of the minimum-mean-square- 
error (MMSE). 

[0033] FIG. 5 is a flow diagram of the method for timing estimation. 

[0034] FIG. 6 is a flow diagram of the method for channel estimation. 

[0035] FIG. 7 is a flow diagram of the method for timing estimation, channel estimation 
and cancellation. 

[0036] FIG. 8 is a block diagram of the present invention. 

[0037] FIG. 9 is diagram of the timing estimation results of the present invention. 

[0038] FIG. 10 is a diagram of the channel estimation results of the present invention. 

[0039] FIG. 1 1 is a diagram of a typical channel impulse response. 

Detailed Description 

[0040] As described in detail in, " Modeling and Simulation of Discretized Data 

Transmission in Very High-Speed Digital Subscriber Line", Dogu Arifler, Ming Ding, 
and Zukang Shen published in EE382C-9 Embedded Software Systems, Spring 2002, 
DMT is a form of multi-carrier modulation. In DMT modulation, available bandwidth is 
divided into multiple narrowband channels. These channels can be viewed as 
frequency-indexed sub-carriers that are modulated and demodulated independently. 
To eliminate inter-symbol interference in DMT, it is desirable to have the bandwidths 
of the sub-channels sufficiently narrow. Substantial research has shown that such a 
system can provide very high bandwidth efficiency with proper allocation of bits and 
energies to these sub-channels. An efficient bit allocation algorithm would 
appropriately assign more bits to the sub-channels with higher signal-to-noise ratios 
(SNR). 

[0041] 

A typical DMT system utilizes a transmitter inverse Fast Fourier Transform (IFFT) 
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and a receiver forward Fast Fourier Transform (FFT). Ideally, a frequency domain 
equalizer following the receiver FFT can correct the channel frequency distortion. 
However, the delay spread of the channel in the beginning of the receiver FFT block 
contains inter-symbol interference from the previous block. As this interference is 
independent of the current block of data, just the frequency domain equalizer cannot 
cancel it. The typical solution cyclically pads a block of prefix data in front of the FFT 
data block on the transmitter side before the block of FFT data is sent to the D/A. This 
block of prefix data is referred to as the cyclic prefix. The prefix data is a replica of 
the last few samples of the FFT data block. On the receiver side, the received signal is 
windowed to remove the cyclic prefix data. If the length of the channel impulse 
response is shorter than the prefix length, intersymbol interference from the previous 
FFT data block is completely eliminated. Frequency domain equalizer techniques are 
then applied to remove intra-symbol interference among DMT subchannels. However, 
since the channel impulse response varies and is dependent upon each individual 
channel, there is no guarantee that the length of the impulse response is shorter than 
the prefix length. Using an adaptive time domain equalizer often shortens the length 
of the channel impulse response. 

[0042] 

Referring now to FIG. 1 , which shows the system level block diagram of a DMT- 
VDSL modem. In practice DMT exploits the power of Discrete Fourier Transform (DFD, 
and it is efficiently implemented using Fast Fourier Transform (FFT). In FIG. 1 , a typical 
modem system of the present invention is described. The input bits are received 1. 
During modem start-up, certain carriers are used as pilot tones 5. Prior to 
transmission both transmitter and receiver know the data sent over these tones. Pilot 
tones are helpful for timing information extraction as well as for channel estimation. 
The bits are then modulated to provide N complex valued symbols 10, typically 
Quadrature Amplitude Modulated (QAM) symbols X( k), 0 <, k <> /V-l modulate N 
orthogonal sub-carriers use IFFT 15 and generate N samples of time-domain signal x 

( n) as in FIG. 2, 2-1. In DMT, a set of QAM symbols X™ are generated as described 
in FIG. 2, 2-2. Due to the conjugate symmetric property of the input data, the discrete 

output signal x will have real-valued samples when the Inverse Fast Fourier 
n 

transform (IFFT) 15 is performed at the transmitter. Prior to the output signal's 
transmission into the channel, a cyclic prefix is inserted at the beginning of the real 
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DMT time symbol 20\o combat Inter-Symbol Interference (IS I) between adjacent 
symbols. If the length of the cyclic prefix is larger than the length of the equivalent 
discrete time channel or the channel impulse response (CIR), then ISI can be effectively 
eliminated. However, in practical channels the CIR can be equivalent to a large number 
of samples, when this is the case, synchronization is severely delayed. To overcome 
large CIR situations, a receiver equalization scheme comprising a T-tap time domain 
equalizer (TEQ) 80 to shorten the length of the CIR and a frequency domain equalizer 
(FEQ) / 10 consisting of A/1-tap filters are utilized to correct for signal phase rotation 
and signal amplitude attenuation. To achieve channel partitioning and to avoid inter- 
symbol interference (ISI), DMT/OFDM extends the modulated transmit vector to satisfy 
the condition in FIG. 2, 2-3. The result is that the transmit vector is extended by 
copying the last ^samples of the multicarrier symbol at the beginning of the DMT 
symbol. This is called cyclic prefix (CP) 20. Before transmitting, the signal is up- 
sampled 30to match the sampling rate of the D/A converter 25 at the analog front- 
end. The converter transforms discrete-time samples of the signal into an analog 
signal, which is transmitted over the communication medium. The transmission 
medium, telephone line channel, introduces noise 35 40 in addition to attenuation, 
distortion and phase rotation to the signal 35 40. The received signal with the 
Channel Impulse Response (CIR) denoted by vector h is shown in FIG. 2, 2-4. This 
denotes the effects of noise 40 and attenuation due to channel effects 35 . The signal 
arriving at the receiver is down-sampled 55 at the A/D converter 45 . The cyclic prefix 
is removed 85 and the signal is demodulated using FFT 90. The resulting QAM 
symbol is as shown in FIG. 2, 2-5. Therefore, the overall system equation in 
frequency-domain can be written as shown in FIG. 2, 2-6. 



Since transmitter and receiver do not share the same clock signal, in the 
beginning they are not synchronized. In order to decode the information properly, the 
receiver must initially synchronize itself with transmitter. Receiver synchronization is 
generally performed in two steps timing offset estimation 50and timing offset 
correction, through symbol clock correction 75 and delay-rotor property 100 . 
Essentially, synchronization aligns FFT windows of both at the transmitter and the 
receiverThe present invention utilizes Maximum Likelihood (ML) estimation to extract 
timing offset information. Since it is a non-data aided (NDA) method, the receiving 
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side can perform timing estimation before FFT, and therefore synchronization can be 
achieved faster. 

[0044] Referring now to FIG. 5 and FIG. 1 , according to a method of the present invention 
a DMT frame of length N+LP is transmitted 125 . Let A/be the number of samples in a 
DMT frame, and LPthe length of CP. The frame of N+LP samples is received at the 
receiver 130 . Since the CP is generated using last LP samples of total N samples, it 
yields maximum correlation when aligned with those last LP samples of the DMT 
frame. This property is utilized to estimate the timing offset information at the 
receiver 140 . At the receiver, the correlation function is calculated for each value of 
estimated delay as in FIG. 3, 3-1. However, in the presence of severe interference 
and/or low Signal to Noise Ratio (SNR), a large number of frames are required to 
extract optimal timing information. The final result is a longer start-up phase for the 
modem. Therefore, rather than using just LP samples for the calculation of Maximum 
Likelihood (ML), the present invention uses the whole received DMT frame of N+LP 
samples to correlate with a pre-stored frame of the same size at the receiver. With the 
same accuracy, this approach helps to shorten duration of the timing estimation, and 
therefore reduces start-up overhead. Further, prior to correlation, interpolation is 
performed on the received frame to estimate integer as well as fractional timing 
offsets 135 . The integer estimate is referred to as the symbol timing offset ( A ) and 
the fractional estimate is referred to as the sample timing offset ( e). Averaging the 
( 4 + e) over multiple DMT frames during start-up provides a very accurate estimate 
of the timing offsets. 

[0045] Fig. 9 shows the performance of the ML algorithm with the use of pre-stored 
frames of length N+LP as taught by the present invention. Simulation parameters 
were, total 2048 sub-carriers, 256 samples CP, VDSL channel as shown in Fig. 1 1 , 
D/A sampling frequency 35.328 MHz. Therefore, the data coming from the DMT 
modulator is interpolated from 1 7.408 Msps to 35.328 Msps using an up-sampler and 
a 22-tap half-band filter. For simplicity, all the sub-carriers were assumed carrying 
64QAM symbols. At the receive side output of A/D is decimated to 1 7.408 Msps. Fig. 
9 also shows comparison of actual estimate with optimum values. 

[0046] 

After the receiver estimates the timing offsets, it should advance or delay its 
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sampling clock by the same amount to get synchronized with the transmitter. This 
procedure is called timing offset correction. The receiver does correction for symbol 
timing offset ( A ) 65\r\ time-domain through delaying or advancing the symbol clock 
75 by integer number of samples 145 . And compensation for the fractional sample 
timing offset ( e) 70 is done in frequency domain 145 through the delay-rotor 
property 100. 

[0047] The receiver also performs channel estimation so that it can compensate for the 
effects of the telephone line channel. Prior art systems are adapted to perform 
channel estimation after synchronization. The present invention provides a novel 
solution that allows for the simultaneous estimation of timing offsets and the channel 
response. The present invention model uses certain carriers as pilot tones, and the 
transmitter sends known QAM symbols over these pilot tones 125 . At the receiver 
130 this information is used to estimate the effect of the channel. The present 
invention utilizes one of the optimal sets of pilot tones. Within a total of A/subcarriers, 
the optimum set of L pilot tones are selected as shown in Fig. 3-2. 

[0048] Assuming the channel to be time-invariant over a DMT frame and the noise to be 
white gaussian, overall input-output relationship of the equation shown in Fig. 2-6 
can be written as shown in Fig. 3-3. Using the minimum mean square error (MMSE) 
criterion the estimated channel impulse response is given 760 as shown in Fig. 3-4. 
The normalized mean square error (MSE) in channel estimation is defined as in Fig. 4- 
1. 

[0049] The comparison of the actual channel with its estimation is shown in Fig. 10. MSE 

-4 

was found to be 2.07 x 1 0 or approximately -36 dB, which is quite 
normalized ... 
low, and hence the estimation is very accurate. 

[0050] estimated channel impulse response has a total of L samples. This process of 

channel estimation is repeated over a burst of DMT frames, and the impulse response 
is found by averaging them. It helps to reduce the effect of AWGN, and reduces MSE. 
After the channel estimation is completed, the modem synthesizes a 1 -tap frequency 
domain equalizer (1 FEQ) 170 , one complex valued tap per carrier. It can be 
implemented by multiplying the FFT output with a complex multiplier (tap). Essentially 
in absence of noise, it is a zero forcing (ZF) equalizer, and the individual tap values 
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are given by Fig. 4-2. And when AWGN is present, according to MMSE criterion the 
individual tap weights can be calculated as shown in Fig. 4-3. 

[0051] Although 1 FEQ is straightforward to implement, it can efficiently compensate 
most for channel attenuation and phase rotation. It should be noted that a circular 
convolution problem is encountered with 1FEQ. This can be easily understood by 
referring to the respective operations of the channel and the equalizer. The channel 
effects are a linear convolution of the incoming signal and the CIR. In contrast, the 
1 FEQ attempts to invert the channel effects after the FFT. It is therefore a circular 
convolution operation. In response to the circular convolution problem, the channel 
impulse response is padded with ( N-L ) zeros 165 , which makes the CIR period the 
same as that of the incoming DMT frame. 

[0052] Fig. 8 shows an illustrative example of a system embodied by the present 

invention. A transmitter 175 transmits the DMT frames as previously described 180 . 
The receiver receives the frames 185 . At the receiving side, an interpolator 195 and a 
correlator 200 provide the timing offset estimations to the symbol clock correction 
circuit 215 and the rotor delay correction circuit 210 . The received frames are acted 
upon by the channel impulse response estimator 790 and the estimation is fed to the 
synthesizer 205 to enable equalization. 

[0053] As described, the present invention discloses a new useful and nonobvious 
solution for DMT modem initialization. The present invention discloses a unique 
combination of timing offset estimation and channel estimation. This unique 
combination provides an improved system and method for fast link initialization 
between modems without significantly increasing complexity in the system 
architecture. 

[0054] (t wj || t h us b e seen t h a t the objects set forth above, and those made apparent 

from the foregoing description, are efficiently attained. Since certain changes may be 
made in the foregoing construction without departing from the scope of the invention, 
it is intended that all matters contained in the foregoing construction or shown in the 
accompanying drawings shall be interpreted as illustrative and not in a limiting sense. 
It is also to be understood that the following claims are intended to cover all of the 
generic and specific features of the invention herein described, and all statements of 
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the scope of the invention which, as a matter of language, might be said to fall there 
between. 

[0055] Now that the invention has been described, 
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